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I. INTRODUCTION 



Determination of the eigenspectrum of the Helmhohz equation for different boundary conditions and geometries 
has been a long standing problem both in physics and in engineering. A particular class of the problem relates to 
finding out the eigenvalues when the corresponding eigenfunctions satisfy the Dirichlet or the Neumann condition 
on the boundary. In fact, one can also have mixed type of boundary conditions, although it comes up less 
frequently in practical situations than the previous two cases. The problem has been pursued in different contexts 
in physics and in engineering. Hence, the importance of solving the Helmholtz equation for different boundary 
conditions can hardly be overemphasized. The standard wave equation that one encounters once too often in 
different realms of theoretical and applied physics reduces to the Helmholtz equation after one separates out the 
temporal part. For instance, the canonical example where one deals with the wave equation is the vibration of 
membranes and plates. The problem crops up often even in acoustic and in mechanical engineering. Perhaps 
the best known example of the Neumann boundary condition is the analysis of the propagation of the TE modes 
in a waveguide. Another common example (in 1 dimension) is that of heat conduction through a rod when 
energy is supplied at one of its ends at a constant rate. Examples of the Dirichlet boundary condition include 
the propagation of the TM modes of electromagnetic waves within a waveguide and the vibration of membranes. 
However, in all these cases simple analytic solutions are available only for a very restricted class of boundary 
geometry. In general obtaining a solution to the problem for an arbitrary boundary can be a formidable task. 
However, many physical situations require that we solve the equation on a general domain. In particular, the 
design of waveguides with a cross-section perturbed from a square may turn out to be useful in eliminating the 
losses due to the corners. The same problem presents itself, although in a different guise, in quantum mechanics 
where one solves the Schrodinger equation to find out the energy states of a particle confined by an infinite 
potential well in 2 dimensions and enforces the Dirichlet condition on the boundary of the well. A variant of the 
same problem is the study of quantum billiards which, however, involves more subtleties. Another area which 
has witnessed a flurry of activities of late is the study of quantum dots. The dots are usually taken to have a 
circular symmetry. But in practice that can hardly be guaranteed. There are bound to be small departures from 
exact circular symmetry. In such a scenario a very natural extension is to consider the confining region to be 
a supercircle [l[ and investigate the resulting spectrum. Recently, Bera et al (H have proposed a perturbative 
approach to the problem where the correction terms were given by a power series expansion and the method 
was applied for the supercircular boundaries. This then calls for a programme to solve the Helmholtz equation 
for a general boundary. Chakraborty et al Q suggested a general recipe which was in the same spirit of Bera 
but furnished the perturbative corrections in a closed form for the case of the Dirichlet condition. On a more 
optimistic note, we can even attempt to proceed the other way round and expect that a knowledge of a part of 
the spectrum may offer an window to probe the exact shape of the dots. In fact, attempts have already been 
made exploring these issues The question 'Can one hear the shape of a drum?' was first posed by Kac as 

an inverse eigenvalue problem to the Helmholtz equation ff\. In the context of quantum mechanics this implies 
whether a knowledge of the energy states of a particle confined in a 2 dimensional infinite well is enough to 
shed light upon the boundary of the well. A counter-example was constructed by Gordon et al [8] who gave a 
pair of isospectral but non-congruent domains. However, by exploiting the inside- outside duality of quantum 
billiard problems [9l-[l3j it has been shown recently that the shape of isospectral billiards can indeed be inferred 
by studying the exterior Neumann scattering problems 14]. Thus, considering its immense applicability, it is not 
surprising that the problem of solving the Helmholtz equation attracts the attention of physicists even today. 

However, the problem, both in its classical and quantum incarnations, is amenable to exact analytic treatment 
only in some special cases. Such problems are generally tackled by invoking the method of separation of variables 
in a suitably chosen coordinate system. Thus rectangular and circular boundaries can easily be handled. The case 
of the triangular boundary is more involved [Tsj . The problem of elliptical boundary, though solvable in principle, 
is quite non trivial. This is because in elliptical coordinates the separation of variables leads to the Mathieu 
equation whose solutions are the Mathieu functions. Mathieu functions have a 2-parameter dependence which 
makes them quite complicated to handle. For any other sort of boundary we virtually run out of a suitable choice 
of coordinate system. The problem becomes even more untraceable for boundaries having no simple geometrical 
shape. This is a major handicap for an analytical study of the 'irregular boundary' problem. An alternative 
way then is to estimate the eigenvalues by numerical means. In fact, till now, most of the efforts at finding out 
the eigenvalues of the Helmholtz equation for an irregular boundary have been along this direction [T6l - [3l| . The 
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analytical aspect, however, has not received much attention except for some occasional attempts. The^ problem 
has been taken up in different contexts by Rayleigh j32j] , Fetter and Walecka 13311 , Morse and Feshbach [34] , Parker 
and Mote [s^] , Nayfeh , Read [s^ and recently by Wu and Shivakumar [sl] and Molinari (39i] . Dubertrand et 
al have employed a similar scheme for studying the propagation of electromagnetic waves in open dielectric 
systems. More recently, Amore '^T] has obtained a systematic approximation to the ground state energy and 
the wavefunction on a general domain in a nonperturbative way. The general spectrum has also been obtained 
employing standard perturbative tools. However, most of these papers address only the Dirichlet condition. In 
this paper we extend the method of [3] to encompass the case of the Neumann condition. This is crucial since 
it gives us an analytic handle to estimate the eigenspectrum in an effective and systematic way. Besides obvious 
academic interest this is also necessary if one wishes to consider TE mode propagation of waveguides through 
channels of arbitrary cross-section. The validity of the method is confirmed by comparing the analytical results 
with those obtained numerically for supercircular and elliptical boundaries. Further, we also propose a semi- 
empirical formula, involving only a single parameter, giving the eigenvalues of the Helmholtz equation for the 
case of the supercircular boundary with both type of boundary conditions. To our knowledge, this is the first 
instance where a semi-empirical formula has been put to use to determine the eigenspectrum. This equips us with 
a simple yet effective way to find out the eigenspectrum for the supercircle. The validity of the formula is judged 
by estimating the energy levels for several low lying modes and comparing them with numerical calculations. 

The paper is organised as follows: In section 2 we set up our general scheme and in section 3 we apply it to the 
case of a supercircle and an ellipse. In section 4 we describe the supercircle and the curious approximate duality 
with respect to the supercircular exponent. The details of the semi-empirical formula for both the Dirichlet and 
the Neumann boundary conditions are provided in section 5. A short conclusion and a discussion of the results 
are presented in section 6. 



II. PERTURBATION ABOUT THE EQUIVALENT CIRCLE 



The homogeneous Helmholtz equation on a 2 dimensional fiat simply connected surface S reads. 



(1) 



with the Dirichlet condition ij) 



on the boundary dT, or the Neumann condition — 



^ - on as, where ^ 

an ' an 

denotes the derivative along the normal direction to 9E. We are interested in finding out k'^{= u) for the case 
of the Neumann condition. Once we are handed any periodic function r{9) = r{9 + 27r) defining the boundary 
of the domain S in flat 2D where we intend to solve the equation we first construct a circle of radius i?o which 
respects the 'equal area' constraint 



As = 



1 



'{B)de = ttrI 



(2) 



where As is the area of the domain S. The next task then is to Fourier expand r(6) about the 'equal area' radius 
Rq at different orders of smallness (denoted by A), 



ri9) = Ro 

where the Fourier series at order a is given as 



(7=1 



oo 

/(-)(0) = ^(C(-) cos 710 



Urn 



sin?if 



n=0 



(3) 



(4) 



This essentially implies that we are seeking a perturbative solution to the problem in the form of corrections to 
the solution for a circular boundary. While the choice of the smallness parameter A is certainly dictated to a large 
extent by the nature of the boundary we are considering, a certain amount of discretion is also required on our 
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part to ensure the efficiency of the perturbative scheme as discussed in detail in the Conclusion of 3] . Of course 
the Fourier coefficients appearing in ([4]) are not all free but constrained by the 'equal area' condition, 



oo (T— i 



(5) 



n=0 v=l 



In particular we have, 



OO 



)2 



(6) 
(7) 



Following Rayleigh [32| we argue that at the zeroth order level the eigenvalue loq of the system confined by r{9) 
will be given by that corresponding to the 'equal area' circle of radius Rq 



ujQ = —J such that p[ j = k^jRo (8) 

Rq 

with p'l ,■ being the j*'' node of the derivative of the V"^ order Bessel function, i.e. J/(pJ j ) = 0. To fine tune the 
expression for the eigenvalues the next step will be to incorporate the correction terms arising due to the deviation 
of the given boundary from the equivalent circle. We do a perturbative expansion of the eigenfunction ^ and the 
eigenvalue uj in terms of the parameter A as 

V' = -00 + AV'i + AV2 + (9a) 



LO = bjQ + AcJl + A W2 + 



(9b) 



where ■00 is the unperturbed eigenfunction with corresponding eigenvalue wo. Plugging in these expansions in ([T]) 
and collecting separately the coefficients of different powers of A produces the following set of equations 



0(0) 
0{\) 



(10a) 
(10b) 
(10c) 



The correction to the eigenfunction at each order will be constrained by the Neumann boundary condition 

VV' • n = (11) 

where n is the normal at the boundary given by r = r{9). Taylor expanding the expression in (|lip about r = Rq 
and making use of (|9ap yields the constraints at each order, 



0(0) 
0(A) 

0(A2) 



^o(i?o) = 0, (12a) 
i^[{Ro) + i?o/^'Vo(^o) - -^F^oiRo) = 0, (12b) 

V'^(i?o) + i?o/^^V'i'(^o) + i?o/^'V^,'(i?o) + ii?^/(i'Vo"'(i?o) - -^P^MRo) 

I Ho 

- -^P^oiRo) + ■^f^''>F^o{Ro) = 0. (12c) 
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In the above expressions prime (/) and hat (') denote differentiation with respect to r and 9 respectively. In the 
following we discuss separately the cases I = and 1^0. 



A. Non-degenerate states {I — 0) 



For the I = state we have, 

^0 = NJoip) (13) 

where Jo is the 0*'' order Bessel function, and iV is a constant which can be evaluated by normalising the 
eigenfunction suitably, ujq is obtained from ([8]) with I — 0, and an appropriate j. 

To correct the wavefunction upto first order we need to solve (jlObp with the restriction imposed by (jl2bp . The 
most general solution to (jlObp is 

oo 

^pi = y^(ap cospO + ap smp6)Jp + uqJq — f^—^NJi (14) 

where the last term is the particular integral. Constraining this solution by demanding that it obeys (jl2b|l and 
matching the coefhcients of the sine and the cosine terms separately provides the coefficients ap and hp 

ap = P^,,iVC(i) Jo(Po,,)/j;(Po,,), (15) 
ap-p;,,^-iV4i)jo(p^,,)/4(/^Oj), P^O (16) 
wi = 0. (17) 

Simultaneously, we also note that the first correction to the eigenvalue vanishes for the Z = case. So any possible 
correction to the eigenvalue can only stem from the second or higher order calculations. We have one more 
constant ag to determine which can be fixed by normalising the corrected eigenfunction over the region S. But 
we will not evaluate it in the present case since the exact value of oq will not be necessary for the calculation of 
the eigenvalues. 

The second order calculation will mimic that of the first order. Now we need to solve (llOcp with wi set equal to 
0. The most general solution is 

oo 

V'2 = V(&p cospO + hp siiip0)Jp + boJo - ^NJi (18) 

where again the last term represents the particular integral. Imposing the constraint (ll2cp on the second order 
correction provides expressions for the coefficients bp and bp together with the second order correction to the 
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eigenvalue, 



2 =^£[c«- 



/ Mp'oj) 
fo,j- 



(19a) 



Np'o.jMPoj) 



+ 



r^(l) ^(1) , c.(1) , /^(l) I c(l) c(l) c.(l) c.(l) 

/ . ["-"p+fc^fc ~^ ^p+k'^k + '-'Ifc-pl'-'fe + '^fc-p'^fe '^p-k'^k 



fe=l 



(Po.,) 



(19b) 



'^piPo,j/ 



\ oo 

'Poj-^oiPojJ V- [^,(1) ^(1) ^(1) o(i) , ^(1) cj(i) c(i) r^(i) , 

4J'(^p'^ .) Z^l'^P+k^k ^p+k'^k '^^\k-p\^k '^k-p^k '^p-k^k 



+ci'' {ip'i, + kM;ik - (po:, - - s%)}] (i9c) 

So at the second order we expect a nonvanishing correction to the eigenvalue in general. The form of 0^2 which 
we have obtained is completely general and holds for any type of boimdary which departs from a circular one by 
a small amount. So for a particular boundary one only needs to calculate the Fourier coefficients Cn^ and Sn^ 
to compute the second order correction to the eigenvalue. As before, the remaining constant 60 can be fixed by 
normalising the wavefunction upto the order of A^. 



B. Degenerate states (Z =^ 0) 

The I ^ states come in 2 varieties, 

For these states we shall assume that = for all a. This is done only to simplify the calculations and will 
not mask any crucial physics. Also for definiteness we prefer to work with one variety, namely 

i^o =NiJi{p) COS W & ""0 = ^ ' '^'i(pU)=^ (21) 
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and just give the result for the other case. The first correction to the wavefunction is obtained as a solution to 
()10b|) and respecting (|12bp 



■01 = apJpC0sp6 + ( aiJi --NiJi+i ] cosW. 

^ — ^ \ LJq L ' 



(22) 



As in the case of the I = case here also the constraint will give the expressions for Op and alongwith the first 
correction to the eigenvalue 



Wo 



2 4ip'u)' 



NiMp'i.j) 



iPLj+pi)ci'ii + ipi:j-pi)ci;Li\ 



(1) 



for p^OJ 



(23a) 

(23b) 
(23c) 



with ai being obtained from the normalisation condition. However, we will not evaluate it here. Unlike in the 
I = case now we expect a non-zero correction to the eigenvalue at the first order itself. The second order 
corrections can be found out in exactly the same manner, using (|10cp & (|12cp . 



h T „ T 

Op Jp CLp 'Jp-\- 1 

^0 ^ 



cos p9 



^2 
UJq' 



pNi 



cosW, 



(24a) 



Pl, + I' 



^21 ^ ^^0 



(2) 



Pi, 



Pi 



3 / n=l 

\ oo 



-2(1 



^2 



) , ^(1) 



2 1 



2p, . — 

'^1-3 n=l 



JLiPu) 



(24b) 



The constants bm can also be determined as mentioned earlier. The equation for ^2 (|10cp also involves ipi in one 
of the inhomogeneous terms. So it may seem that we need to know the coefficient a; (which we have not found 
out here) to determine the form of ip2- However, it turns out that the terms involving ai cancel out amongst 
themselves so that we really do not need to explicitly evaluate it. For the other variety. 



^Po = NiMp) sinW 



(25) 
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we can proceed in an analogous manner. However, for the sake of completeness we give only the final results here, 

ui p 



and 



Wo 
a„ : 



\^ apJpSinpO + ( a;J; --NiJi+i ] sin/0. 



= a 



(1) (pm+1 



21 



(1) 



p=l 



h T - - T 

^P'-^P n^P P~^^ 

Wo Z 



for p 7^ I 



W2 



— -rJi+2 Ji+i 

wg 4 Wo 



sin Z6', 



(26a) 

(26b) 

(26c) 

(26d) 



W2 _ 1 w2 / - 2f 



, 1 /pf,-3/2\ - 



4 \ - P 
^ \ Pi,j ' / „=i 



(1) 



+21 ^|2;-ri|^ 



(r./ + 2/2)Ci^_/,,-(n/-2/X!- 



2/| 



2 



OO 



) ^(1) 



71=1 



2 1 



11=1 



2p; 



^ OO 



'j' n=l 



2 Jn(p;,j) 



(26e) 



III. APPLICATION TO SIMPLE CASES 



A. Introduction 



In the preceding section we have formulated the general formalism in detail. In this section we put this formalism 
to test by determining the spectrum when the boundary is supercircular or elliptical. Direct comparison with 
numerical results are also made. We have solved the equation numerically, using the Partial Difi^erential Equation 
Toolbox"^^ of MATLAB® by defining our supercircular or elliptical boundary into it. 



B. Supercircular boundary 



Supercircle is a special case of a superellipse [4l| whose equation is given by the Lame equation. 

It L,|t 



1. 



(27) 
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with t > and rational, a and b being positive real numbers. In the literature they are also known as Lame 
curves or Lame ovals |42l |. Superellipses can be parametrically described as, 



X =acos^/*(u), 
y =5sin^/*(u). 



(28) 
(29) 



A supercircle will obviously correspond to setting a — b. For each value of 'i' - the supercircular exponent we 
obtain a different curve. The shapes of the supercircles for different values of t are shown in FIG[T] Evidently, 
t = 2 corresponds to a circle whereas t — 1 corresponds to a square with its sides rotated by an angle of j. The 
case t ^ oo also corresponds to a square with its sides parallel to the axes. In these limiting cases the Helmholtz 
equation can be trivially solved for both Dirichlet and Neumann type boundary conditions without invoking any 
special techniques. We will, however, be interested in a general value of t. More specifically, in the spirit of the 
perturbative formalism we have developed we would like to find the spectrum when t slightly deviates from 2. 
The problem was addressed in [3| for the Dirichlet condition. In this paper our aim will be to solve the same 
problem for the Neumann type boundary condition. When t > 2 we shall take only the real positive values of 
cos* (u) and sin~(w). Moreover, while calculating the eigenfunction we will confine our attention to the region 
^ u ^ ^ and exploit the symmetry of the region to continue to the other quadrants. In polar coordinates the 
equation for the supercircle is, 

(30) 



(cos* 6* + sin* 6')i/* 



and the radius of the equal area circle is. 



Rq = a 



[r(f)] 



(31) 



We define the deformation parameter to be A = < — 2. The Fourier expansion of r{9) is given as 



r^Ro 



1 + A ^ C[l^ cos 4ne + X'^J2 ^in cos AnO 



n=0 



(32) 



with 



C^^' is given by equation ^ to be -0.0017552 and C, 



^in — 



1 



AniAn"^ - 1) ' 



(2) 
4 



32 



= 0.0357983 [|. 



A knowledge of these few coefficients is sufficient to determine the eigenvalue corrections to second order. In 
FIGlU FIG [3] and FIG|3]we have plotted the first few levels (continuous line) using our scheme and compared 
with the numerical results (discrete points). We have varied the deformation parameter in the range — 1 ^ A ^ 1. 
We find that the analytical results tally with the numerical ones for small perturbations around the circle, i.e. for 
small value of |A|. Within the range 1.5 < t < 2.5 the matching is quite impressive. As we increase the deviation 
the agreement becomes poor which is not unexpected. Even then for a distortion as large as t = 3 the agreement 
with numerical results for many of the low lying modes is quite encouraging indeed. 



C. Elliptical boundary 

As a second example we determine the eigenvalues when the boundary is an ellipse. The case of elliptical 
boundary has been studied extensively in the literature because it is one of the few cases where the Helmholtz 
equation can be solved exactly. Even then it is an arduous task since it involves the Mathieu functions and 
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very often one has to depend upon numerical methods j23l-l25l|. Recently, Wu and Shivakumar [38| attempted 
an analytical solution to the problem. Here we will use our perturbative algorithm to solve the problem for the 
Neumann type boundary condition. The polar equation for an ellipse with semiaxes a and b is 

= . (33) 



l-(l-^)cos^ 



We will choose, 



a — b , ^ 

(34) 



as our deformation parameter although the eccentricity might seem to be a natural choice. This is one instance 
where one needs to exercise some discretion to ensure the effeciency of the perturbation theory and enhance the 
precision of the analytical predictions (see the Conclusion in Q). (|33p can be recast as, 

r = Ra[l + A cos 261 - ^A^ + ^A^ cos 461 + 0{X^) + ...]. (35) 

A comparison with our general Fourier series yields the Fourier coefficients relevant to our purpose, C2^'' — 1, 
c(^)=-iandCf = |. 

The results for the elliptical boundary are shown in FIG H] and FIG|6l The deformation parameter varies from 
-0.25 to +0.25. Here also the analytical predictions tally with the numerical calculations for a fairly wide range of 
A. The only notable exceptions are the cases where the levels repel each other. For those cases also the formalism 
works well for one of the levels (J; sinW) involved. However for the other level (J; cosW) it falters. For the other 
cases, the matching between analytical and numerical results is quite pronounced between —0.15 < A < 0.15. 
For a larger deformation the agreement becomes less satisfactory. It should also be mentioned that in some 
of the higher states our numerical estimates may not be very accurate. For example it is seen that while the 
analytical results forecast exact degeneracy of some of the higher states throughout the entire range of A, the 
numerical technique that we have used reflects that degeneracy only in a certain range of A. Outside that range 
the degeneracy is lifted. So the apparent disparity between the analytical and the numerical estimates has got 
to do with the limitation of the numerical scheme we have employed rather than the failure of the perturbative 
formulation for the higher modes. Perhaps one can employ some other numerical techniques and see whether the 
agreement improves for the higher states and whether the degeneracy is preserved. 



IV. THE SUPERCIRCULAR DUALITY 



In this section we discuss a very interesting duality for the case of a supercircle which will prove to be immensely 
helpful when we construct the semi-empirical formula in the next section. Without loss of generality henceforth 
we choose a = 1 in the equation of a supercircle. We also restrict t > 1, as before, so that we always have a 
convex curve. The area of the supercircle is given by. 

We show that for indices t and t' related by the following 0, 

11 

- + - = 1: 1 < < < 2 and c» > > 2, (37) 
t t' ' - - - - : \ ! 



the shapes of the supercircles are almost identical. Clearly, a circle is self dual since for a circle t ^ t' = 2. 
On the other hand, the 2 squares — one corresponding to i = 1 and the other corresponding to i' = cxd — are 
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also dual to each other. But the areas of the 2 squares are different and differ just by a scale factor. However, 
the duality we define in this way is not exact. In fact, the two aforesaid cases are the only ones where we have 
an exact duality and our naive argument of scaling the area (the scale factor being unity for the self-dual case) 
holds true. Nevertheless, it is an enticing prospect to explore how far this simple scaling argument can give us 
an estimate in the other cases. For our purpose it is convenient to recast (using the properties of the Gamma 
function, T{2z) = 2l!_r(z)r(z + 1/2)) the area A of the supercircle (with 1 < t < 2) as 

_ 4V¥r(i + i) 

" 2VT(i + i)- 

For 2 < i' < oo the area when expressed in terms of its dual exponent t defined by the relation (|37p is. 

In FIGlHwe have shown the major axes, AC and PR, for t — 1.5 and its dual parameter t' — 3. The corresponding 
minor axes are QS and AC respectively. In general, the ratio, a, of the major axis to the minor axis (for example 
AC/QS) for 1 < t < 2 is 2^t~2\ The ratio, a', of the major axis to the minor axis bv 2 < t' < oo (for example 
PR/AC) is 2(^-f). Now it is easy to see that the above two ratios are same, i.e. a = a' , ii t and t' satisfy the 
duality relation (|37p . a also measures the scale (the ratio of two major (or minor) axes, i.e. PR/AC (= AC/QS)) 
through which the supercircle is inflated under the duality transformation (|37p . In case of exact duality the area 

should scale by a^, i.e. A{t') — a^A{t). In FIGHwe have plotted A = ° ^A(t')^''* against t to see whether such 
a scaling law gives a good estimate of the area for the dual figure. The maximum aberration from the scaling 
law is seen to take place at two points viz. tmin ~ 1.2485 and tmax ~ 5.0242 and it is obvious that tmin and 
tmax satisfy the duality relation (j37p . The deviation is zero at t = 1 and t = 2 thereby vindicating our previous 
assertion that these 2 cases exhibit exact duality. For the other cases the deviation is very small (within 3.3%) 
and by virtue of that we can safely exploit the duality relation for extracting the eigenvalues in the domain t > 2 
once the results are known for 1 < t < 2. 



V. THE SEMI-EMPIRICAL FORMULA FOR SUPERCIRCULAR BOUNDARY 



In this section we propose a semi-empirical formula for determining the eigenspectrum for the case of the 
supercircular boundary for both type of boundary conditions. In the following we shall consider the Dirichlet 
and the Neumann boundary conditions separately. The case of the Dirichlet condition is discussed in detail. The 
Neumann condition can be similarly treated and hence only briefly dwelt upon. 



A. Dirichlet boundary condition 



Given a supercircle of area A we consider a square of length L and a circle of radius R and each having an area 
A. The length of each side of the 'equal area' square is, 

L./f^, (40) 

and the radius of the 'equal area' circle is, 



[r(f)] 
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The solution to ([T]), when the boundary is a square (S*) of side L is, 

^ = -sm(-^)sm(^j, (42a) 

and hence the eigenvalue, 

Us = k'.^inl+nl)^ (42b) 

where Ux and rij, are positive integers. It is then easy to see that ip — on dS (i.e. on x — 0^ x — L, y = and 
y = L). When the boundary is a circle (C) of radius R, the solution is. 




with pij = kcR and 



V'=| V2 , ^cos/g\ |, (43a) 



2 

= = ^ (43b) 



where pij is the j*'' zero of the Z*'' order Bessel function, J; (p). Now to build up a semi-empirical formula for the 
case of a supercircle we will attempt to find an expression for the eigenvalue w(t) in the range 1 < t < 2. Further, 
the duality discussed in the preceding section will allow us to extend our results to 2 < t < oo making the whole 
domain as 1 < t < oo. It is easy to figure out from FIG[T]that in the range 1 < t < 2 the boundary has a shape 
which is an intermediate between a circle and a square. So it makes sense to try an expression of the form, 

uit) - f{t)LJs + (1 - /(i))c^c, (44) 

i.e. we take an average of Wg and luc with suitable weight factors. We note that all the eigenvalues are calculated 
for boundaries enclosing an equal area given by (j36p . Hence the L and the R in the equation (|42[) and (|43l) are 
identified with the ones appearing in the equation (|40l) and (j4ip respectively. To find an explicit form of f{t) we 
note that in the limiting cases t = 1 and t = 2 the semi-empirical formula must reduce to, 

a;(l) = LUs and w(2) ~ ujc- (45) 

In other words, in such limits, we must have, 

/(1) = 1 and /(2) = 0. (46) 

Clearly, there can be infinitely many choices of f{t) fulfilling the above criteria (|46p . The simplest choice is 
f{t) = 2 — t. Unfortunately the eigenvalues calculated with this choice does not match very well with the 
corresponding numerical results. So an exact form of these weight factors cannot be found out from the above 
arguments only and is best obtained through a 'trial and error' method. We have explored different choices of 
f{t), calculating the eigenvalues and comparing them with the numerical results. As the next plausible candidate 
we take f{t) = (2 — t)" and try to estimate the best n using numerical results. Finally we have settled for the 
choice, n — 3, so that 



f{t)^{2-tr. 



(47) 
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So the final expression for the eigenvalue for a supercircle with an exponent t (1 < t < 2) is furnished through 
the following expression; 

u;it) - [(2 - tf£s{n,,ny) + [1 - (2 - tf]£cipu)] (^8) 

where £s{~ L'^lOs) and £c(= ttR'^lOc) are the eigenvalues for a unit area square and a unit area circle respectively. 
For a particular excited mode nx,ny and pij have to be chosen properly. In Table U] we have shown the exact 
mapping between the parameters Ux and Uy of £s and pij of £c for the first 10 states. The eigenvalues (including 
the degeneracies) £s and £c are arranged in increasing order of magnitude and are then matched one to one. Note 
that all the states with I ^ for the circle are doubly degenerate and hence are written twice. Each pair of and 
Uy corresponds to one pi j. Here we have excluded the trivial case Ux = Uy — 0; Ux = 0,ny = 1; = l,ny = and 
pi i — 0. Having done this we next embark upon the task of finding the eigenvalues in the regime 2 < i' < oo. To 
find an explicit form of f'(t') in this region we note that in the limiting cases t' = 2 and t' = oo the semi-empirical 
formula should yield, 

U!{2) — and uj{oo)—uJs, (49) 

i.e. we require, 

f'(2) = and /'(oo) = 1. (50) 

This is where the duality relation (|37|) comes into play. Plugging it in ()47|) we arrive at the following functional 
form of fit'), 

Clearly, such a choice of f'{t') satisfies the limiting conditions. So, finally for t' lying in the range 2 < i' < oo the 
duality implies that the expression for the eigenvalues will be: 



a;(t') = 



1^ 3 
t 



1-F 



£s{nx,ny) + il - [- }£c{pL, 




mi) 



(52) 



Using the above expressions we have evaluated the eigenvalues for the first 21 states (including the degenerate 
ones) for a supercircle and compared them with the numerically obtained values. The comparisons are shown in 
FIGEI FIG|9]and FIGfTOl We find that the agreement is excellent over a wide range of the supercircular exponent 
t except for a few cases which we discuss in the next section. 



B. Neumann boundary condition 



In case of the Neumann boundary condition we must have = on 9E. The vanishing of the normal 
derivative implies that for the square boundary with side L the eigenfunction should now be, 

2 fnxTTX\ fnyTTy\ 

This, however, has no effect on the eigenvalues with the exception of Ux — 0, riy = 1 and = 1^ Uy — for which 
one now has non trivial tp- On the other hand, for the circular boundary of radius R the eigenfunction remains 
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the same as before, but the Neumann boundary condition dictates that we must have, 

-kl = P^ (54) 

where pj ^ is now the j*'* zero of the derivative of the l*^ order Bessel function, i.e. J'lip'i j) — 0. Barring these 
shght modifications, the basic recipe for finding out the eigenvalues remains unaltered from the Dirichlet case. 
The eigenvalues are given by the expressions ()48p and ()52p with properly matched Ss and £c- In Table ITTl we 
have illustrated the exact matching between the parameters by considering some low lying modes. Here also 
we have excluded the trivial CcLSG Tlx — T^y = and /9g_i = = as it was done for the Dirichlet case. We 
have plotted the eigenvalues for the first 18 states (including the degenerate ones) and compared them with the 
numerical results in FIGlTT]and FIG112I Once again the matching is found to be reasonably good save and except 
a few levels whose discussion we postpone to the next section. Comparison of the empirical and the numerical 
eigenvalues for the first seven states for a supercircle with exponent i = 1.5 for both type of boundary conditions 
is given in Table IIIII 



VI. CONCLUSION 

In this paper we have generalised the formalism developed in Q to encompass the case of Neumann type 
boundary condition. We have also explicitly verified the results for two boundary geometries - supercircle and 
ellipse. This is quite significant since the Helmholtz equation is not exactly solvable for such geometries. Here we 
briefly highlight the salient features which make our formalism different from the existing ones. The scheme we 
develop is completely different from those prevalent to attack such kind of boundary value problems. However, it 
does have a parallel in the Rayleigh-Schrodinger type potential perturbation in quantum mechanics. In fact, one 
of the primary motivations of the order by order construction of the Fourier series in 3] was to make contact with 
the time independent perturbation theory in quantum mechanics. In 2-level systems in quantum mechanics one 
witnesses the phenomena of level crossing when the Hamiltonian is perturbed. This is most easily achieved by 
placing the system in a constant electric or magnetic field 43] . The analogy with quantum mechanics is further 
consolidated with the observation that boundary perturbation can also induce crossing of the levels. Although 
in quantum mechanics one usually solves a Dirichlet boundary value problem, here we show that the phenomena 
of level crossing can appear even in Neumann type problems. Another unique feature is that unlike any other 
perturbation scheme ours does not require a separate formalism to treat degeneracies. Degenerate states are 
handled equally elegantly as the non-degenerate ones. As opposed to earlier schemes where the perturbation 
parameter is an artificial one here a real one, extracted from the parameters defining the irregular boundary, 
serves as a perturbation parameter. It is found that for any perturbation the higher states are more affected. We 
have plotted beyond the 20th excited level and have shown that the analytical perturbation procedure developed 
here does reasonably well even for such higher states. Finally, our scheme is a one shot process. If there is an 
equation for the family of closed curves (like that in an ellipse or a supercircle) one can have the full spectrum 
just from a knowledge of the parameters of the equation. We exploit the parameters of curve equation as the 
perturbation parameter which is a measure of the deviation from a circle and obtain the results in one go. One 
has to just substitute the parameter values to get the new spectrum whereas in all other existing schemes one 
has to find each time the Fourier coefficients and then substitute even if the curve is of the same family. 

Besides the phenomenon of level crossing, another interesting case is when the levels just manage to evade a 
crossing and veers off after touching each other tangentially. For the 2-level quantum mechanical system mentioned 
here one can in fact set up a criteria to predict whether or not a pair of levels will cross. It will be interesting to 
explore the possibility of formulating an analogous criteria for boundary perturbation. A closely related problem 
is to investigate why the perturbation theory falters whenever the levels display repulsion for elliptical boundary. 
We have seen that for the elliptical boundary, the agreement between the analytical and numerical results is quite 
outstanding for small values of A. The only blemish is the inability of the scheme to predict the level repulsions. 
A proper understanding of this level repulsion can be a potential topic for further study. 

We have also put forward a single parameter semi-empirical formula giving the eigenspectrum of the 2 dimen- 
sional homogeneous Helmholtz equation with Dirichlet or Neumann condition on a supercircular boundary. It 
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is found that the prescription gives results to a high degree of accuracy for a wide range of the parameter t for 
both type of boundary conditions barring a few cases. In fact we have checked expUcitly for 1 < t < 5 and the 
agreement is rather impressive throughout the range. This is quite remarkable considering the fact that a semi- 
empirical formula involving just a single parameter, the supercircular exponent t, yields such accurate results. 
This assumes more significance because the problem is not exactly solvable analytically. Even the perturbative 
scheme that we have developed works only for small deformations about the circle. On the other hand, the semi 
empirical formula, supplemented with the duality relation gives us an analytic handle to extract the eigenvalues 
even when t ^ 2. A possible generalisation would be to attempt a similar formula when t < 1, i.e. we have 
a concave curve. We would also like to point out a very peculiar feature that has plagued the semi-empirical 
results: the abnormal behaviour of the Bessel function J2- The semi-empirical formula for the Dirichlet boundary 
condition does not provide the anticipated results whenever there is a contribution to lOc from the zeros of J2 
which is quite mysterious indeed. The first such case arises from /52,i which contributes to the eigenvalues of the 
second and the third modes. It is clear from Table IIIII that for the Dirichlet case the agreement between the 
empirical and the numerical values is outstanding and the maximum error is 0.6% except for the second and the 
third states where the error is within 4%. The next victims are the ninth and the tenth modes which receive 
contribution from the second zero p2,2 of J2 and in both cases the degree of mismatch is little pronounced but the 
relative errors are within 5%. In the case of the Neumann condition also there is disparity between the numerical 
and the semiempirical predictions for some of the eigenvalues. The first instance occurs in the case of the first 
and the second modes both of which involve a contribution from P2 i . Table IIIII also shows that the errors for the 
aforesaid modes are around 9% for t — 1.5. For the other modes the matching is very good and the errors are 
within 2%. The next levels afflicted with this anomaly are the eleventh and the twelfth ones which involve P2 2- 
This mysterious behaviour of J2 is something which cannot be accounted for at this point and calls for deeper 
study. Strangely enough, all these J2 states also involve level repulsion. In this context, it may be mentioned that 
Chakraborty et al 0] had earlier reported the failure of perturbative treatment in the cases where the levels repel 
each other for the case of Dirichlet condition on an elliptical boundary and there also the levels were associated 
with J2. Here too the analytical scheme for the Neumann boundary on the supercircle reveals that the J2 states 
display level repulsion. From the results at hand it seems as if level repulsion and J2 abnormality go hand in 
hand. But we do not know whether this is sheer coincidence or there is an underlying connection between the 
two. In any case we think it would be worthwhile to look into this aspect. 
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Table I. Matching of the parameters of £s and £c for finding u){t) for the Dirichlet boundary condition 
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Table II. Matching of the parameters of £s and £c for finding a;(t) for the Neumann boundary condition 
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Table III. Comparison between empirically and numerically calculated eigenvalues for t = 1.5 
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Figure 1. Shape of the supercircle for different values of t. 
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Figure 2. Comparison of the eigenvalues obtained numerically and analytically for a supercircular boundary (in units of 
■^) with Neumann condition for the first 11 states. 
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Figure 3. Comparison of the eigenvalues obtained numerically and analytically for a supercircular boundary (in units of 
■^) with Neumann condition for the states 12 to 20. 
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Figure 4. Comparison of the eigenvalues obtained numerically and analytically for a supercircular boundary (in units of 
with Neumann condition for the states 21 to 26. 




Figure 5. Comparison of the eigenvalues obtained numerically and analytically for an elliptical boundary (in units of 
with Neumann condition for the first 7 states. 
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Figure 6. Comparison of the eigenvalues obtained numerically and analytically for an elliptical boundary (in units of 
with Neumann condition for the states 8 to 15. 
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Figure 7. Justification of the duality relation. 
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Figure 8. Comparison of the eigenvalues obtained numerically and empirically for a supercircular boundary with Dirichlet 
condition for the first 10 states. 
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Figure 9. Comparison of the eigenvalues obtained numerically and empirically for a supercircular boundary with Dirichlet 
condition for the states 11 to 15. 




Figure 10. Comparison of the eigenvalues obtained numerically and empirically for a supercircular boundary with Dirichlet 
condition for the states 16 to 21. 
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Figure 11. Comparison of the eigenvalues obtained numerically and empirically for a supercircular boundary with Neumann 
condition for the first 9 states. 




Supercircular exponent (t) 

Figure 12. Comparison of the eigenvalues obtained numerically and empirically for a supercircular boundary with Neumann 
condition for the states 10 to 18. 



